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ABSTRACT 

Little Higgs models have an enlarged global symmetry which makes the Higgs 
boson a pseudo-Goldstone boson. This symmetry typically contains sponta- 
neously broken U(l) subgroups which provide light electroweak-singlet pseu- 
doscalars. Unless such particles are absorbed as the longitudinal component 
of Z' states, they appear as pseudoscalars in the physical spectrum at the 
electroweak scale. We outline their significant impact on Little Higgs phe- 
nomenology and analyze a few possible signatures at the LHC and other 
future colliders in detail. In particular, their presence significantly affects 
the physics of the new heavy quark states predicted in Little Higgs models, 
and inclusive production at LHC may yield impressive diphoton resonances. 
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1 Introduction 



Little Higgs models [1] have recently emerged as an alternative solution to the naturalness 
problem of the Higgs sector in the Standard Model (SM). In these models, the origin of the 
electroweak (EW) scale is identified as new dynamics in the multi-TeV range that involves 
the spontaneous breaking of some extra continuous symmetry. While the exact nature of this 
dynamics, the UV completion, is undetermined — scenarios involving strong interactions [2], 
iterated Little Higgs models [3], or supersymmetry [4] have been proposed — the low-energy 
effective theory below the UV scale, A, is determined by the assumed symmetries. Various 
realizations of the Little Higgs symmetry structure have been proposed [5,6,7,8,9]. 

If at the scale A an exact global (i.e., ungauged) symmetry is spontaneously broken, the 
spectrum contains Nambu-Goldstone bosons (NGBs) which are massless and have no renormal- 
izable interactions. In analogy to low-energy QCD, we denote the would-be decay constant of 
these scalars by F. Naive dimensional analysis relates this scale to the cutoff A by F ~ A/4?t. 

In Little Higgs models, the SM Higgs doublet is among these NGBs. This would be a natural 
explanation for a weakly interacting Higgs sector. However, since the Higgs doublet does have 
nontrivial renormalizable interactions — gauge, Yukawa, and self-couplings — there must be 
interactions in the initial Lagrangian which break the global symmetry explicitly. In general, 
such symmetry-breaking interactions introduce a one-loop Coleman- Weinberg potential [10] 
and a Higgs mass proportional to m 2 H ~ A 2 /(167r 2 ) ~ F 2 . Without fine-tuning the parameters 
or adding additional fields, one then derives F ~ v, i.e., the new symmetry-breaking scale is 
near the electroweak symmetry breaking (EWSB) scale [11]. 

The new ingredient in Little Higgs models is the mechanism of collective symmetry break- 
ing, as it was observed in the context of deconstructed extra-dimension models [12]. Each 
renormalizable scalar interaction breaks the postulated global symmetry explicitly, but leaves 
a continuous subgroup intact. Spontaneous breaking of this remaining global symmetry then 
still implies the existence of a massless NGB. However, all interactions together break all global 
symmetries explicitly, and no particle stays massless to all orders. The Higgs doublet is found 
among the scalars that still remain NGBs as long as only a single symmetry-breaking coefficient 
(spurion) is present in the Lagrangian, but acquire masses once all spurions are turned on. For 
instance, for two spurions <7i,<?2 the resulting quadratic term in the effective Higgs potential is 
proportional to (g\ / '4ir) 2 (g 2 / 'An) 2 A 2 . As usual, EWSB is triggered by such a term, and thus we 
have a three-scale model with 

v ~ F/47T ~ A/(4tt) 2 . (1) 

Since the UV-completion scale A is parameterically two orders of magnitude above the EW 
scale, any sign of the associated dynamics is strongly suppressed, and we are left with the 
virtual effects of new particles at the intermediate scale F. For a consistent implementation 
of collective symmetry breaking, enlarged symmetries must be introduced in all sectors of the 
model, so we expect new vector, spinor, and scalar particles with masses of order F. The 
low-energy traces of these particles can be computed and have been used to derive limits on 
the parameter space of any given Little Higgs model [13,14,15,16,17,18,19]. 

In this paper we study the phenomenological consequences of a particular property of the 
Little Higgs mechanism. To allow for an EW doublet among the NGBs, the spontaneous 
breaking of the global symmetry group typically involves a reduction of the group rank, e.g. 
577(3) — > SU(2). While the Higgs doublet in this example corresponds to the off-diagonal 
broken generators (analogous to the kaons in QCD), there is also one broken diagonal generator 
(analogous to the rj). If there were no explicit symmetry-breaking terms, this particle, which 
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acts as a pseudoscalar in its fermionic couplings, would behave as an axion, so we may call it the 
pseudo-axion of a Little Higgs model. Similar particles show up in other models which involve 
enlarged global symmetries, e.g. technicolor models [20,21], see-saw topcolor models [22,23], 
or the NMSSM [24,25]. To determine the detailed model structure, one must experimentally 
establish any pseudo-axions as part of the Higgs sector and determine their properties. We find 
that these states have significant, broad impact on Little Higgs phenomenology. 

2 Pseudo-axions in Little Higgs Models 

In the bosonic sector, two different lines of model-building realize collective symmetry breaking. 
In models built along the lines of theory-space or moose models [1,6,7,8,9], the global-symmetry 
representation is reducible, i.e., in the scalar sector there are several distinct multiplets with 
gauge quantum numbers. The gauge coupling of any multiplet acts as a spurion which reduces 
the global symmetry down to the exactly realized gauge symmetry. This is spontaneously 
broken down to the EW gauge group SU(2) x U(l), and some NGBs are absorbed as the 
longitudinal components of new heavy vector bosons. However, if there are several scalar 
multiplets, there are not enough gauge bosons to absorb all the NGBs. The masses of the 
uneaten linear combinations are proportional to two or more spurions (i.e., gauge couplings) 
and do not appear up to one-loop order. Introducing scalar self-couplings as additional spurions, 
the structure becomes more complicated and more scalar multiplets are needed to keep the Little 
Higgs mechanism working, but the line of reasoning remains unchanged. 

For concrete examples, let us consider the Minimal Moose [6] and Simple Group [7] mod- 
els. In the first, with exactly one scalar coupling turned on, the global symmetry breaking 
is [S77(3)] 4 — > 577(3). The group rank is reduced by k = 6 units, yielding 6 pseudo-axion 
candidates. The rank of the spontaneously broken gauge group SU(3) x SU(2) x U(l) is 
r = 2 units larger than the SM EW group, so 2 axions are eaten to make vectors massive 
and k — r = 4 of them remain in the low-energy spectrum. (Two are EW singlets, while the 
other two are the neutral members of two real EW triplets.) Similarly, in the Simple Group 
model, if exactly one scalar coupling is turned on, the pattern of global symmetry breaking is 
[5f/(4)] 3 -> [SU(3)] 2 x SU(2), while the gauge group is SU(4) x U(l). This yields k = A and 
r = 2, so there are k — r = 2 pseudo-axions in this case. 

As an alternative, some models implement an irreducible representation of the symmetry 
group in the scalar sector [5,8]. In this case, to provide independent spurions in the gauge sector, 
the enlarged gauge symmetry must be such that part of it commutes with the EW gauge group. 
In this setting, for which the Littlest Higgs model [5] serves as an example, scalar self-couplings 
arise at one loop from integrating out the heavy vector bosons and fermions. Here, the global 
symmetry breaking is SU(5) — > 5*0(5), with rank reduction k = 2. If the extra gauge group 
is SU(2), one pseudo-axion (r = 1) is eaten in its symmetry breaking, and one (k — r = 1) 
remains in the spectrum. On the other hand, if the extra gauge group is SU{2) x U(l), as 
originally proposed, both axions disappear from the spectrum. 

This example demonstrates that the new particles become unphysical if all extra broken 
£7(1) symmetries are gauged. However, in that case we expect the corresponding number of 
new Z' bosons with masses of order F. These states are generally easy to detect at future 
colliders, either directly as resonances in qq or e + e~ annihilation [15,26], or indirectly via the 
observation of contact interactions and of mixing with the standard Z boson [13,14]. In some 
models, these effects can be used to rule out much of the parameter space from existing data 
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alone [13,14,15,16,17,18,19]. In the following, we therefore consider the situation where the 
extra groups are ungauged [15,16,17] and the pseudo-axions are physical. 

3 Pseudo-axion Interactions 

Let us consider the case of a single pseudo-axion which we denote by rj. As an EW singlet it 
does not couple to EW gauge bosons at leading order. Furthermore, if rj is an exact NGB, it 
has no potential and does not couple to Higgs bosons either. However, in the fermion sector 
the situation is different. To account for Higgs Yukawa couplings, Little Higgs models contain 
interactions of fermions with the NGB multiplet, inducing Yukawa couplings for the axion. If 
the U(l) symmetry generated by rj is parameterized by 

f = exp —77, (2) 

where F is the symmetry breaking scale of the Little Higgs model, each field ipi transforms like 
V'i ~~ * i^i't) where $ is the corresponding U(l) charge. 

In some models, for some of the axions, these couplings are fixed by the symmetry structure 
that generates EWSB and can thus be deduced from the study of observables in the gauge, 
Higgs, or fermion sectors. In other models they cannot be determined this way. Typically, the 
axion interactions depend on additional parameters, so their observation provides independent 
information on the high-energy structure of the model. In the following, we present three 
specific models where axions play a role and discuss their interactions in some detail. 

3.1 Mass and decays 

For a phenomenological discussion of an rj axion, we need an estimate of its mass. If rj is an 
exact NGB, it would be exactly massless. This appears to be the case in some of the models we 
discuss in the following subsections. However, even then the U (l) v symmetry must be explicitly 
broken at some stage, so that rj picks up a mass; a massless particle coupled to fermions (even 
if this coupling is suppressed by v/F and CKM factors) would induce a long-range force many 
orders of magnitude stronger than gravity and is therefore very strongly ruled out. 

In general, the U(l) v symmetry will be anomalous. Then, if is below the QCD scale, r) 
becomes a Peccei-Quinn axion, which is ruled out for the parameter range of interest. Hence, 
we can assume a lower mass limit in the GeV range. 

In any model, we can also put a rough upper limit on m v by considering its influence 
on the Higgs mass. To this end, we first note that a vertex rfH 2 in the scalar potential 
would re- introduce a quadratic divergence in the Higgs mass at one loop, i.e., a correction 
Am 2 H ~ A 2 /(167r 2 ) ~ F 2 . If there is such a term, its coefficient should be parameterically of 
order l/(167r 2 ), so that the induced Higgs mass correction is at most of order v 2 . Next consider 
rjH scattering: the NGB Lagrangian will contain derivative interactions of the form drfdH 2 /F 2 
which at one-loop order yield an effective rfH 2 vertex with a quadratic divergence proportional 
to (m 2 /F 4 )A 2 /(167r 2 ) ~ m 2 /F 2 . To keep this consistent without fine-tuning, we generally have 
to require m v < v, i.e., the EW scale is an upper limit for the pseudo-axion mass. 

All rj couplings to SM particles are suppressed by v/F, hence the rates for direct production 
are generically a factor of (v/F) 2 smaller than the corresponding Higgs production channels. 
For the decays, throughout the mass range we expect a similar pattern as for a light Higgs boson: 
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dominant (tt,) bb, t + t~ , . . . branching ratios (BRs), depending on kinematic accessibility, and 
some fraction of gauge-boson pairs, i.e., gg, 77, and Z^^y, ZZ^*\ WW^*\ Since the axion is 
CP-odd, the vector-boson pair branching fractions are loop-induced and therefore stay small 
even above the thresholds for on-shell WW and ZZ pair production. 



3.2 The Littlest Higgs model 

Here we derive the pseudo-axion interactions in the Littlest Higgs model [5] with only one 
£7(1) gauged [15,16,17]. While this model has been extensively discussed in the literature, the 
pseudo-axion interactions have so far been ignored. 

In this model, the NGB multiplet is parameterized by a 5 x 5 matrix 

E = (exp ^U] S , (3) 
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and n = — W -A V / V10 h T . (4) 



We have included only the physical states: the Higgs doublet, the 77 singlet which multiplies the 
(canonically normalized) generator diag(l, 1, —4, 1, l)/2\/l0 , and which is a complex triplet 
written as a symmetric 2x2 matrix. The triplet becomes heavy [i.e., its mass is of order F] 
and has little impact on low-energy phenomenology, so we ignore it in the following. 
The kinetic term for the NGB multiplet is given by 

^Tr[(^S)*(^S)] = + i(^) 2 + -.., (5) 

which fixes the 77 field normalization. In writing the Yukawa interaction we have to allow for 
77-dependent factors 

£ = exp , (6) 

where the normalization has been adjusted for later convenience. The third-generation fermions 
are the left-handed quark doublet qi = (iz,,&L) T , the right-handed singlet tn, and the new 
singlets Tr,Tl. We also define the matrices 

(it 2 ^T l iq L 0\ 

Xl=\ -tq T L and %T 2 2 = diag(0, 0, -ir 2 *)/2 . (7) 
\ 0/ 

With these definitions, the Yukawa interaction has the form [5,19] 

C} Y = X.F^HrTt [E*(iT 2 )E* XL ] - X 2 F^f R T L + h.c. (8) 

The parameters A),i,2 are real numbers determined by the differences of fermion 27(1)^ charges. 
The model has no prediction for these numbers (note that anomaly cancellation is not an issue 
since the U(l) v symmetry may well be anomalous), so we leave them as free parameters. In 
particular, we allow for (3 2 7^ 0. If this is the case, the T quark mass is protected by a chiral 
symmetry and is therefore naturally of order F, the U(l) v breaking scale. 1 



1 Previous papers [15,16] assumed that the T quark is vectorlike, leaving this fine-tuning problem unsolved. 
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When the Yukawa interaction in Eq. (8) is expanded in powers of 1/F, it yields a mass term 
for the new fermion T, a mixing of t and T, a Higgs Yukawa interaction of the top quark, and 
higher-order terms: 

C Y = -X.FIrTl - \ 2 FT R T L + ^pl\ 1 t R h T ~q L 

i i (9) 

+ — (2 - A) - A)AiWT L - —p 2 \ 2 T R r]T L + ■■■ + h.c. , 

where q~i = iT 2 qh- 

The mass matrix is diagonalized first by a rotation of the right-handed fields with the 
parameters 



s = Ai/A, c=A 2 /A, A=y / A2 + A|, (10) 

which are expected to be 0(1), and then by subsequent rotations of the left- and right-handed 
fields where the rotation angles are suppressed by powers of v/F. Keeping only the leading 
terms, we write the interactions in terms of the physical fields H and r\ and the quark masses 
M T = \F and m t = sc\v, where v is the Higgs vacuum expectation value: 

£ Y = -M T f R T L - m t t R t L 

+ sc^HT R T L - -—Hf R t L - —Hi R t L n -,x 

t C V V \ ll J 

+ -^(3 V T R T L + i^f3'r)t R T L - ' V i R t L + ••• + h.c. 

CS V V t 

The coefficients of the axion couplings depend on the fermion U(l) v charges as: 

P= -^[(2-p -p 1 )s 2 -p 2 c 2 }, (12a) 

P'= ^Q-Po-fh + fo), (12b) 

P"= ^=[(2-(3 + (3 2 )s 2 + (3 lC 2 } (12c) 

From the mass-eigenstate Lagrangian of Eq. (11) we can read off the vertex structure. By 
construction, rj is CP-odd. While the particular values of the Yukawa coupling coefficients are 
specific to a given model, their orders of magnitude with respect to the mass scales of the 
theory, shown in Table 1, are dictated by the SU(2) transformation properties and thus are 
generic for Little Higgs models. The relation of the TT and tt couplings is a consequence of 
the fact that, after Little Higgs symmetry breaking, the heavy states are vectorlike, while t is 
a chiral fermion. For TT, unbroken EW symmetry forbids a coupling to the Higgs doublet, 
so this coupling must be proportional to v/F which arises from t — T mixing. Similarly, the 
coupling tir) is forbidden for unbroken symmetry since rj is a singlet. The chirality assignments 
of the mixed couplings are model-dependent: If the heavy singlet T is replaced by a doublet, 
Pr and Pl have to be exchanged. For T as a heavy triplet, the structure of TiH is unchanged 
while the Tir] couplings are suppressed by v/F. 

Although not previously defined in the literature, we posit that the Yukawa coupling for 
down-type fermions in the Littlest Higgs model can be written down as: 

C b Y = -^F^bR Tr \£Y iX l\ + h.c. , (13) 
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TTH 0(f;) 

TtH 0(l)V L + 0(f)V R 

UH 0(1) 



TT V 0(1) T5 

Tfy 0(1)^ + 0(f) P L 

fify g(i)75 



Table 1: Chiral structure of the Yukawa couplings in Little Higgs models with a heavy SU(2) 
singlet T. Pr^ are the chirality projectors. For the coefficients, only the order of magnitude is 
indicated. This structure remains valid for the fi model presented in Sec. 3.3. 



where Y\ is one of the hypercharge operators defined in [19], 

y 1 = iliag(3,3,-2,-2,-2). (14) 

The pseudoscalar coupling is then 

g vbh lry b with g vbb = -/jF^A* (15) 

Clearly, (3$ is an additional new free parameter. 

The free parameters in the model are then one Yukawa coupling, Ax, the mass scale F, and 
the four differences of U(l) fermion charge assignments, /? , 1,2,3; A 2 is fixed by the top quark 
mass condition. We consider the case F = 4 TeV, which is motivated by the precision EW 
limits for the Littlest model [19]. 
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Figure 1: Left: Mt (solid) v. Ai in the Littlest Higgs model with F = 4 TeV. Also plotted is 
A2 (dashed), which is constrained by the choice of \\ and the known value m t . Right: Heavy T 
quark branching ratios to Wb (green), tZ/th (blue), andtt] (red), for three choices of the U{1) 1) 
charge differences A), 1,2-' 1,1,1 (solid); 1,0,1 (dashed); 0,0,1 (dotted). 



The presence of the pseudo-axion rj alters the decay spectra of the heavy quark T (see 
Appendix A), drastically for some regions of parameter space. Note in Eq. (11) that the 
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T — > tr] partial width is proportional to fa 2 , which is a pure difference of the fa and contains no 
mixing factors, whereas the T — > tH partial width is protected by 577(2) and therefore appears 
only as a consequence of T — t mixing. The latter is true also of the decays T — > . 
(Ignoring mass effects, T tH = T tz = \T bW .) Due to the sum rule = 1 / Af + I/A2, both 
\i and A2 must be larger than \ t . While there is no general obstacle against \i > X 2 (i.e. 
sin > cos), this is disfavored by large mixings of the third-generation EW currents [19]. So we 
take Ai to be bounded from below by A t , and from above by \/2\ t using the condition Ai = A 2 . 
This translates to a range 0.72 < Ai < 1.02. In Fig. 1 we plot m T v. the allowed values of Ai, 
as well as the consequential A 2 , which closely tracks my. As expected, in the Ai — > Af M limit, 
A2 becomes large and almost totally responsible for tjit, which becomes (perhaps unnaturally, 
although one could argue that Ai = Af M is "natural" as well) much larger than F. For the 
lower extremal limit, t — T mixing vanishes and the partial decay widths to Wb, tZ, tH vanish 
correspondingly (the decoupling limit). The right-hand side of Fig. 1 shows the resulting T BRs 
for three choices of the relevant fa. For large mixing and equal U{l) n charges, BK(tr]) bottoms 
out at around 5% and dominates only when the mixing becomes very small and uit grows quite 
large. For other choices of fa, e.g. only fa 7^ 0, fa triples and BR(t^) dominates everywhere, 
with the T quark total width being O(10 — 100) GeV. For negative /3 ,i, which are in no way 
unnatural, the non-77 decays would be practically squelched for all Yukawa couplings, and the 
total T width can be hundreds of GeV. The shapes of all curves in Fig. 1 are independent of 
the choice of F. 

The dominant decays of the i] are to bb, gg and 77, the latter two being loop-induced 
couplings (see Appendix B), and for m v > 2m t , decays to top quark pairs. While there 
will also be loop-induced decays to ZZ and W + W~, the BRs to observable final states are 
small fractions of the already rare decay rates, so we do not consider these further. There 
are in principle also direct couplings to the lighter fermions, but proportional to the fermion 
mass squared. We ignore these, as the only straightforward, high-efficiency signal would be 
to muons and that BR is typically an order of magnitude smaller than to photons. The BRs 
are independent of the scale F, since it appears the same way in all partial widths; and to 
the percent level, independent of Ai. Table 2 shows the BRs to the dominant final states for 
various values of m^. Note that from m# stability arguments we don't expect m, q > v, but this 
is order-of-magnitude, so we show one case with m v > 2m t for illustration. 



m v [GeV] (fa) 


150 (1,1,1,1) 


150 (1,0,1,1) 


150 (1,1,1,0) 


250 (1,1,1,1) 


400 




0% 


0% 


0% 


0% 


99.3% 


r] — > bb 


38% 


72% 


0% 


17% 


0.02% 


99 


61% 


29% 


99.7% 


83% 


0.5% 


77 -> 77 


0.17% 


0.08% 


0.27% 


0.27% 


0.002% 



Table 2: Dominant branching ratios of the r\ in the Littlest Higgs model. The case fa = (0, 0, 1, 1) 
is identical to (1,0,1,1). For the casem v = 400 GeV, the Pi assignments are irrelevant, as decay 
to top quark pairs overwhelmingly dominates. 

It should be obvious that the phenomenology of the Littlest Higgs model has an extreme 
dependence on the presence of the pseudo-axion and the pattern of U(l) v fermion charge 
assignments fa, which are undefined. It would not be inaccurate to say that the model is fairly 
incompletely defined. The situation is even more complicated than presented above, since there 
should be some U(l) v charge assignment for every SM fermion, greatly increasing the number of 



7 



free parameters, although one could reasonably argue based on flavor-changing neutral current 
(FCNC) constraints that the charge assignments are identical at least across generations. Even 
the presence of just the third-generation fa introduces enough new parameters to present a 
serious phenomenological challenge at future colliders. The positive aspect of this is that 
measuring these parameters may provide clues to the model's UV completion. Furthermore, 
the Little Higgs mechanism doesn't depend on any of the fa. 



3.3 The fi model 



Recently, Schmaltz has proposed a moose-type model [9] where the Little Higgs mechanism 
is implemented in a very economic way. A fi term which explictly breaks some of the global 
symmetry is responsible both for the absence of fine-tuning and for EWSB. The EW group is 
enlarged to gauged SU(3) x U(l). There are two nonlinear sigma model fields, each of which 
parameterizes a coset space U(3)/U(2): 



$1 = exp 



'■ F *n~ 








1 











$0 



exp 



[->¥_ 

















(16) 



where 2 



e 




F 2 = F 2 + Fl 



(17) 



As usual, the global Little Higgs symmetry is broken radiatively by the gauge and Yukawa 
interactions that trigger EWSB and give the Higgs a mass, but they leave a global U(l) v 
symmetry intact. The Coleman- Weinberg mechanism thus generates a negative Higgs mass- 
squared and a positive quartic coupling, but without U(l) v breaking it does not contribute r/ 2 , 
r] 4 , or rfh 2 terms. These are generated only by the \i term: 



-V = ji 2 $}$ 2 +h.c. = 2F 1 F 2 /j, 2 cos 



( Ft, 
\V2Fu 



2F 2 F 2 v ' 24F?Fi 
Thus, the complete potential up to quartic order is (k = Fi/F 2 + F 2 /F 1 > 2): 

-V = -(5m 2 



(rfh) 2 +... 



(18) 



fl 2 K 2 



12F X F 2 V 4 



if | 3(h)h)rf 



+ ■■■ 
(19) 

Here, 5m and 5X are the one-loop contributions to the Higgs boson mass and quartic coupling 
from the Coleman- Weinberg potential given in [9] . From this we read off the r, mass 



m r 



= s/k/j, > \/2 fi. 



(20) 



To minimize the amount of fine-tuning, \x is chosen to be roughly of the order of the EW scale. 
Also from Eq. (19) we deduce the connection 



m 



2 H = -2(5m 2 + m 2 ) . 



(21) 



2 There are other possible choices for the generator T v that multiplies the r\ field, e.g., T 8 or diag(0, 0, 1). How- 
ever, after EWSB these choices introduce kinetic mixing of the r\ with unphysical Goldstone bosons. Removing 
this mixing by appropriate field redefinitions is equivalent to choosing T v proportional to the unit matrix. 
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The interesting property of this model is that the 77 mass is predicted, and can be calculated 
based on other experimental observables. In principle, only the masses of the Higgs boson and 
either the heavy T quark or one heavy vector are needed to determine the scales Fi and F 2 . 
One can then calculate 5m 2 [9] and predict m v by Eq. (21). However, 5m 2 depends on the 
theory cutoff A and has a noticeable dependence, so in practice performing precision fits to the 
model would be difficult, since the physical impact of the cutoff would not necessarily be clear. 
Nevertheless, within some reasonable uncertainty, measuring these values would be a first test 
of whether or not a discovery satisfied this particular model. 

There are two possible gauge charge assignments for fermions in the \i model, called types I 
and II. The first has heavy partners of the top, charm and up quarks, and all generations carry 
identical gauge quantum numbers. This is, however, anomalous. This model also appears to 
be ruled out from precision E W data [9] , so we do not discuss it further. The second model has 
slightly different charge assignments to eliminate the anomalies, and is not yet ruled out. In this 
scenario the top, strange and down quarks have heavy partners, and the Yukawa interactions 
are given by 

C — - \\t 1: R®\^T,L - \%,R^T,L ~ ^ k *\®p k T , L 

+ --- + h.c. 

where the SU(3) triplet is ~§>t,l — { t2< 1l^l) t and d,u sum over the first two generations; we 
assume the up-quark sector to be generationally diagonal, which is not necessary but simplifies 
the model. We diagonalize the mass matrix first with 0(1) right-handed singlet-field rotations 
h,R — > ct R + sT R , t 2 ,R — > — st R + cT R , with s = XiFi/M T and c = X 2 F 2 /M T . We then rotate 
the left-handed fields by an 0(v/F) angle N x m t jM T (see below). The leading (to 0{v/F)) 
Yukawa terms are 

C= -M T T R T L -m t t R t L 

Ht R t L — g V TT HT R t L + g vtt Ht R T L + Qhtt HT l T r 

- - ( 2 3) 

+ ig-qtt V t R t L + IQHTT V T R t L + 1 77 t R T L + IQrfTT V T R T L 

V 

+ ••• + h.c. 

where m t = \\\ 2 vF '/ '\P1Mt and Mt = a/ Af F 2 + A|F| . The couplings are given by 

m t N 2 m t 2 Aq N imt N 2 m 2 N 3 v 

9r,tt ~V2F~^ : 9r)TT ~ 9HTT ~ + 2M~ T ' (24) 

with the abbreviations 

Nl ~ F* x,x 2 ' ^-"fifT' Na ~ Y 2 • (25) 

Comparing this with the Littlest Higgs model, we observe that here the axion properties 
are fully defined, i.e., no additional parameters have to be introduced. Furthermore, there are 
heavy partners of at least one quark and lepton of each generation. By contrast, in the Littlest 
Higgs model the presence of heavy fermions beyond the top-quark partner T is not necessary. 
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We adopt the parameter set of Ref. [9], consistent with existing EW and flavor data and 
the preference for a light Higgs boson. The complete set of inputs, excluding the lepton sector, 
is Fi j2 , X\ :2 , X b , A M , X d , and /i. To simplify the set we assume A}' 2 <C Xl' 2 = A|, so there 
is essentially no mixing between the SM down-type quarks and their heavy partners, whose 
masses are simply X2F2. Deviations from this matter very little for the gross phenomenology. 
We further choose X\ 2 to minimize m T given F 12 and m t . X b and A" are fixed by m^, m c 
and m u . Our only free parameters are then F 1>2 and ji, with the constraint that F 1 not be 
as small as the EW scale, F > 2 TeV from EW precision constraints (primarily AT and 
four-fermion operators), and F 2 > F\ to avoid too much mixing which would lead to fermion 
non-universality. F\ = T 2 lies right at the edge of the limits on the latter. The cutoffs of the 
theory, Ai ;2 , are nominally 47rFi )2 , but this is somewhat vague. Their precise choice affects the 
predicted Higgs mass, but because of the inherent uncertainty as to their definition we simply 
choose the suggested [9] value of A 1 = A 2 = A = 5 TeV and vary /i over a slightly broader 
range than the calculated mu would suggest is allowed or favored by data. 

Starting then with the "Golden Point" in Ref. [9] of 

F 1 = 0.5 TeV, F 2 = 2 TeV, A = 5 TeV (26) 

we calculate 



m T = 990 GeV, m s , D = 700 GeV, m z , = 1.2 TeV, m w , = 950 GeV . 



(27) 



For ji > 150 GeV, mu is lower than the LEP direct exclusion limit of ~ 114 GeV. For /i < 
120 GeV we find that m H is greater than the SM 95% c.l. upper bound from precision EW 
data. However, in general in Little Higgs models, larger values of m H are allowed because of 
additional positive AT contributions from the other new content. We therefore do not restrict 
ourselves to a lower bound on fi. For this Golden Point, as fi — > 0, mu plateaus at around 
450 GeV. To illustrate the dependence on A, if instead we explicitly set A = 4nFi then only 
ji > 135 GeV is allowed. It is noteworthy that m v < 2m t for this parameter choice. 
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Figure 2: Left: m v v. /i for four choices of the ratio Fi/F 2 : 1 (solid), 1/2 (dashed), 1/4 (dot- 
dashed), and 1/6 (dotted). Middle: mu v. ji for fixed F 2 = 2.0 TeV and various choices of 
F1/F2: 1 (solid), 1/2 (dashed) and 1/4 (dot-dashed). Right: mu v. fi for fixed F\f F2 = 1/4 
and various F\ [TeV]: 0.5 (solid), 1.0 (dashed), 1.5 (dot-dashed), and 2.0 (dotted). 
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Figure 3: 77 branching ratios in the [i model for the Golden Point, as discussed in the text, as 
a function of m v . 



Moving away from the Golden Point, we observe several interesting features. First, while 
m v depends only on the ratio F\jF 2 rather than their individual magnitudes, and grows larger 
for smaller ratios, for smaller ratios it is quite easy to violate the LEP direct exclusion limit 
on m H . The trend for increasing ratio is shown in Fig. 2, left panel. The upper end of each 
curve is cut off by the condition that mn > 114 GeV. In general as F±/F 2 approaches about 
1/3, m v > 2m t becomes possible for /i > 200 GeV. 

The right panel shows the (more complicated) trend for mn- As F\ moves closer to F 2 , 
m H prefers to become quite heavy, although as both values become very large, this trend can 
reverse and uih becomes lighter. For example, m# plateaus at around 800 GeV for small 
fi, although at the extremal point Fj = 3.3 GeV the range of mu suddenly drops radically 
to around 100-200 GeV. For Fi > 3.4 TeV, EWSB does not occur, providing a strong upper 
bound on the allowed values of Fi. For smaller ratios Fi/F 2 , both the low-/i plateau for mn 
and the upper limit on fi become smaller. The latter is cut off by the Higgs becoming lighter 
than the LEP direct exclusion limit. A more thorough understanding of the allowed parameter 
space based on EW precision data would require a complete EW fit, which is beyond the scope 
of this work. We would still expect the EW data to prefer a light Higgs, but the upper limit 
could be considerably larger than in the SM, which is why we do not impose an upper limit on 
m H in Fig. 2. 

Decays of the 77 are more straightforward in this model than that of Sec 3.2, because all 
couplings are precisely defined. This is true for the bbq coupling as well, which in this model 
we derive to be 

which is exactly like the tir] coupling (Eq. (23)), except that there is no term proportional to 
iVi from mixing since there is no heavy partner of the b quark. It has the interesting properties 
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that it appears only at higher orders in the expansion of $1,2 (yet remains 0(v/F)), and that 
it is driven rapidly to zero as Fi and F 2 approach each other. Couplings to muons and tau 
leptons will have the same form. The coupling to charm quarks will technically have the same 
form as that for top quarks, but here the mixing term is trivially small due to the ratio m c /Mc- 
We can then straightforwardly predict the partial widths for rj to decay into SM particles. 

However, there is one complication: unlike in the Littlest model, here there is a ZHrj 
coupling which arises from the kinetic term for $ 1; $ 2 (note that both Higgs triplets have the 
same SU(3) W x U(l) x quantum numbers). 

C = (D^D^ + (D^D^ 2 (29) 

with the covariant derivative expressed by the physical fields: 

ZeA. + gfil-tDZ, 0\ 
D^ k = d^ k + i[ $, + ... fc = l,2. (30) 

\ 0/ 

The dots stand for the part of the covariant derivative containing the five remaining heavy gauge 
bosons, which we leave out here. Plugging in the expansion of the nonlinear sigma-model Higgs 
fields, we derive a ZHrj coupling 

CzHv = ^ N2Z ^ H - Hd ^) ( 31 ) 

with N 2 defined previously. Such a coupling is forbidden for unbroken EW symmetry (77 is a 
singlet) so it must be proportional to v/F, but this is partially compensated for by large N 2 
for the Golden Point. 

We plot the BRs for 77 to decay into SM particles, including the ZH final state, in Fig. 3, as 
a function of for the Golden Point. Decays to the final state bb dominate for m v < 200 GeV. 
Above this gg dominates until the ZH threshold is crossed at around m v = 290 GeV. Note 
that BR(77) remains non-trivial over an extremely large range of /i. We can make a rough 
comparison with how a typical two-Higgs doublet model pseudoscalar behaves at low values 
of tan/?. Here, the role of tan/5 is played by the ratio of $ decay constants, Fi/F 2 . Due to 
constraint from lepton universality, we must have F\ < F 2 , or tan/3 < 1, which is the opposite 
of what is required in e.g. the MSSM. 

Although we do not plot the decay BRs for the Higgs, it is worth pointing out that for 
much of the parameter space where the 77 does not decay to ZH, instead the decay H — > Zr\ is 
kinematically allowed and constitutes a very large partial width. We discuss the impact of this 
on general Higgs sector phenomenology further in Sec. 4. 

3.4 The original Simple Group model 

The Simple Group model introduced in Ref. [7] displays an axion phenomenology that is a 
synthesis of the two previously discussed models. The low-energy effective theory contains 
two Higgs doublets. The global symmetry structure exhibits breaking of the semisimple group 
[S , f/(4)] 4 — > [.St/^)] 4 , while one SU(4) subgroup is gauged and broken down to the SM gauge 
group SU(2). Neglecting heavy singlet scalars, the four fundamental scalar multiplets can be 
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approximately written as 



e +iH / F e +i ^/ F (0,0,F,0) T , 
e- iH / F e- iE "/ F (0,0,F,0f, 
w e +iH / F e +i ^/ F (0,0,0,F) T , 
M> 2 ^ e- iif / F e- iEd / F (0,0,0,F) r . 



(32) 
(33) 
(34) 
(35) 



The matrices are 
/0 



V2 



Q g hd 



hi 







E c/ = ^diag(l, 1,-3,1), £ D = ^cliag(l,l,l,-3). (36) 



V hi 0/ 
We introduce the abbreviations 



£ u/d = exp[ir] u/d /(6F)} 



to obtain 



($ 2 )« 
(* 2 )< 



r 3 



£ 3 



1 

72' 



^»,3 H T^hujdij 



^ 3 ^y^h-ujOij 



l 



-F5j,4 H T^hdjSij 



(37) 

(38) 
(39) 
(40) 
(41) 



Considering the potential necessary to break EW symmetry, we can construct four bilinears 
$1^1, $2^2) ^ > l x I / 2 5 and $2^i- I 11 those with like indices the Higgs doublets cancel, so these 
terms approximately vanish. The potential for r] is generated by the remaining terms 



V v = 2b l2 Re 



+ 26 2 iRe 



(42) 



namely, 

$1*2 ~ F 2 (Udf (exp[-2itf/F]) 34 = -^ d f(hlhdU , $2*i » ($ t i*2) t • (43) 
The coefficients 6^ are expected to be of order F 2 . The expansion in powers of 1/F yields 



V v = 2(612 + &21) 



Re [/i^/id] - 



% + % Im KM _(!^ R e KM 



2F 



(44) 



Note that the r] couplings prefactor cannot vanish if EW symmetry is to be broken. Introducing 
the standard Higgs- field components (cf. also [16]), we get (as usual in 2-Higgs-doublet models, 
the scalar Higgs bosons are denoted by h and H, the pseudoscalar by A) 



Re [h{h d ] = - [( Vl + Hc a - hs a ){v 2 + Hs a + hc a ) + c^s^A 2 + 2H + H~)] (45) 

(46) 



2 
A 



Im [h\h d ] = — [vxcp - v 2 s p + Hc a+l3 - hs a _p] 
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We observe that a potential is generated only for the linear combination r/ + = (r) u + r]d) / a/2. 
This pseudo-axion is analogous to the one in the \i model. However, the mass of i] + is very low, 
of order v 2 /F. A further effect of the above potential is mixing of rj + with the pseudoscalar A 
of the second Higgs doublet, suppressed by v/F. 

The orthogonal combination r/_ = (r) u — r)d) / y/2 remains massless, similar to the 77 of the Lit- 
tlest Higgs model. Without introducing extra symmetry-breaking terms we have no prediction 
for m n _ . 

Due to the mixing interaction, there is an hrj + r] + vertex with a coupling of the order v 3 /F 2 . 
This modifies the standard Higgs width by an amount of 



1 / Ami v 5 15 

r "-^~mf-lt^~IFm iMeV< (47) 

and gives a BR of order 5 — 10% into an r] + pair. This decay may be detectable when measuring 
Higgs BRs at a future linear collider. 

The quark Yukawa Lagrangian with r/'s present is 

£y = ^iXi,rXi,lF& + *2X2,rXi,lFZI + hx 3 ,RXi,LF£ 3 d 

+ AiV^Xi,^)^ - A 2 v / 2i(x 2 >i< ? L)^ - X 2 V2i(x 3 ^ d q L )e d • (48) 

Ignoring the coupling Ai (light quark limit) as well as A^ (for the down-type quarks) [16], and 
inserting the vacuum expectation values of the Higgs fields and keeping at most trilinear terms, 
we obtain the fermion mass matrix 

/ o \ / t L \ 

(Xi,r,X2,r,X3,r) A 2 ucos/3cos 7 ^ \ 2 F& • (49) 

\A 3 usin/3sin^d \ 3 F$J \x2,lJ 

Here the additional angles 7 and <fi parameterize different scales within different Higgs mul- 
tiplets. One linear combination of X.2,r an d Xz,r ls ^ ne right-handed top quark, while the 
orthogonal linear combination mixes with Xi,l and X2,l to give two heavy quarks with masses 
of the order F. Here, the physically relevant rotation is between the left-handed fermions. 
The resulting structure of the rj couplings is determined by the general rules given in Sec. 3.2. 
Regarding the couplings, the r] + is something like a pseudoscalar Higgs rendered ultralight 
by a sort-of Little see-saw mechanism between the scales v and F. In contrast, the is a 
pseudo-axion in the pure sense, with similar properties as described in the previous section. 

The extreme complexity of this model is manifest. We therefore do not try to examine its 
phenomenology in detail. Rather, we wish to point out generically that the addition of more 
heavy quark and lepton partners can multiply the partial widths of the 77 to gg and 77, which 
would in general increase the gg — > rj — > 77 rate at hadron colliders, discussed in the next 
section. 



4 Phenomenology 

We concluded in the previous section that the 77 coupling structure to SM fermions and gauge 
bosons resembles the couplings of a CP-odd Higgs boson. Compared to standard two-Higgs 
doublet model interactions, however, there are three important differences: (i) all couplings to 
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SM particles are suppressed by a common factor v/F, which reduces direct production rates 
by (v/F) 2 , although for gg — > r\ at hadron colliders this may be compensated by the presence 
of multiple heavy particles running in the loop; (ii) the Ttr/ coupling is not suppressed and can 
be similar in magnitude to Tth, altering the phenomenology of the new heavy quarks from that 
expected by strict SU (2) symmetry; (hi) the ZHrj interaction is not allowed in all models. 

4.1 Production related to heavy T quarks at hadron colliders 

The pseudo-axion rj has essentially Yukawa couplings to the SM fermions, which are large 
for the top quark. Thus, one would anticipate that associated ttr] production would occur, 
analogous to top-pion [27] or two-Higgs doublet model pseudoscalar production [28] but the 
(v/F) 2 suppression factor renders this channel useless by trivial comparison to the SM Higgs. 
While the rate for TT is not suppressed by (v/F) 2 , it is unfortunately hugely phase-space 
suppressed: the TT rate is already sub-fb level. 

However, the T quark can decay to tr\, often with significant BR and sometimes dominantly. 
At a minimum the T quark phenomenology is changed in a non-trivial way, since the i] does 
not participate in SU(2). Because the coupling is highly model-dependent, and not even 
completely defined in some models, predictions rapidly become exercises in having too many 
free parameters. On the other hand, precisely measuring all BRs, both to the "standard" final 
states tH, tZ, and bW [17], and to trj, could provide deep insight into the full structure of the 
model that is realized in nature, if the Little Higgs mechanism is discovered. 

Heavy T quarks can be produced at hadron colliders singly or in pairs, analogous to SM t 
quark production. However, in general the single-T rate is larger than that for TT because of 
the enormous phase-space suppression in heavy pair production. Consequently, only single-T 
phenomenology has been studied seriously [29]. This study found that even observing TeV scale 
T quarks in a clean channel is difficult with large luminosity, and the final state tH may be 
impossible except at the luminosity-upgraded SLHC [30]. While the single-T production rate is 
different in the \x model compared to the Littlest Higgs model that Ref. [29] investigated, it is 
similar. Because of the extremely complicated final states requiring multiple large, complicated 
backgrounds, we do not attempt to study any of these potential signals in detail here, deferring 
this to a later publication [31]. We do, however, outline some general features one can expect 
in the various models. 

For the Littlest Higgs case, EW precision data already limit F > 4 TeV, in which case 
single-T production rates are fewer than 10 events per 300 fb _1 . This is far too few events to 
utilize, although at the upgraded SLHC [30] it might be feasible to observe the T quark in Wb 
decays [29]. The same study showed that the decay tH was barely feasible for tjit = 1 TeV, so 
by quick comparison of rates we predict it would be hopeless to see either the tH or trj modes of 
a Littlest Higgs even at SLHC. Unfortunately, while the case Ai ~ A t for Littlest Higgs makes 
BK(ti]) dominant, for this parameter choice Mr is simultaneously driven well above the scale 
F: for F = 4 TeV, this region of Ai yields My ~ 20 — 40 TeV, so that the TT final state is not 
even accessible at LHC. Estimates of the signal cross section alone suggest it would not even 
be accessible at a 200 TeV VLHC. However, it is worth investigating the non-extreme case of 
Ai ~ A2 for a VLHC, which is more likely in any case from the point of view of fine-tuning 
of the model and the T quark's role in cancelling the Higgs mass quadratic divergence coming 
from the SM t contribution. 

The situation is different for i] in the fi model, however. For a type II model, EW precision 
data do not constrain the scale of new physics all that much: T 12 = 0.5, 2 TeV is in good 
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agreement with current data. For the allowed range 11 < 150 GeV, ttih ranges from 450 down 
to 120 GeV, m v from 10 up to 310 GeV, BR(T -> tr/) is roughly 1/3, and BR(r/ -> 66) is 
from 90% to 27%. The heavy T quark mass remains fixed (here, 990 GeV) so the rate to a 
triggerable final state is often large. The main issues will be that m v is large enough to avoid 
the large QCD continuum 66 contribution, and that it also does not overlap the Higgs signal. 
Due to finite m b i resolution, this means that the resonances will have to be separated by about 
50 GeV to be separately visible. 

In the \i model there is also the interesting prospect that r\ can decay to ZH, if its mass 
is large enough, typically > 300 GeV. This would allow for the extremely unique decay 
signature T — > trj — > tZH, which is admittedly an extremely complicated final state with low 
efficiency for identification, but also has much smaller backgrounds than any of the standard 
T quark final states, or T — > trj — > tbb. For an rj mass above the tt threshold ZH remains the 
dominant decay mode, with a branching ratio of three quarters, while tt is nearly the remaining 
quarter, so the ZH signal likely remains feasible for large m v . 

4.2 Direct production at hadron colliders 

Another possible production mechanism, which appears to be more interesting than that from 
T decays, is gluon fusion, which proceeds via the axial 27(1)^ anomaly, analogous to the rf 
anomaly of U(1)a in QCD. In general, the more heavy fermions in the model, the more the 
suppression factor v/F is overcome in the loop-induced coupling to gluon pairs, and thus the 
higher the production rate. As in SM Higgs phenomenology, the only viable decay channels at 
the LHC will be weak boson final states, of which 77, similar to gg — > H — > 77 [32,33], offers 
the best prospect due to high efficiency from not having subsequent decays. 

We calculate the rates of gg — > rj production in the Littlest Higgs and /i models for a 
range of allowed parameter choices. Our results are shown in Fig. 4, where we plot both the 
NLO signal absolute cross sections and differentially bin the continuum diphoton background, 
which includes the direct, 1- and 2-fragmentation contributions at NLO [34] 3 . We have applied 
kinematic cuts of Pt(i) > 40 GeV and |?7 7 | < 2.5, as well as an efficiency factor of e 7 = 0.8 for 
each photon to be identified in the detector. 

For the Littlest Higgs model, we select the minimum allowed scale choice, F = 4 TeV, four 
choices of /3j which reflect various ggr] and bbrj coupling strengths, and we consider the mass 
range 100 < m ri < 350 GeV. Unfortunately, these signals appear to be invisible everywhere. 
The most optimistic point is for the case $ = [1, 1, 1, 1] just below the tt threshold, where at 
the SLHC with 3000 fb _1 the signal might be about 4a. Above 350 GeV, the rj would decay 
for all practical purposes only to tt final states, which would be lost in the tt continuum that is 
more than two orders of magnitude larger. The other choices of /3j shown exhibit qulitatively 
different, but rather moot, behavior. For = 0, the cross sections are typically an order of 
magnitude smaller than for the choice f3\ t 2 — 1, as f3\ directly affects the tir] coupling. As m v 
increases, the ratio X1/X2 becomes more important. As expected, the cross sections are larger 
for p 3 = compared to non-zero values, as the partial width to 66 goes to zero, allowing for a 
larger BR(77). Negative values for /3 ,i, which result in dominant T — > trj decays as discussed 
previously, do not result in an enhanced signal here. The reason is that the tir] coupling remains 
approximately the same size, but turns negative, resulting in a cancellation between the t and T 
loops. The choice $ = —2, —2, 0, 1 produces signal cross sections approximately the same size 

3 Diphox NLO continuum background distributions courtesy of Thomas Binoth. 
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Figure 4: Diphoton signal for gg — > r] — > 77 at the LHC, shown as individual points which 
reflect the total signal cross section for each parameter space choice. The symbols for the 
Littlest Higgs points are shown in the plot legend along with the chosen Pi sets. The symbols for 
the n model points are shown with their corresponding choices for Fi,F 2 . Two data points are 
shown for each value ofm v in the Littlest Higgs model, representing the extremal allowed values 
of X±: At and \f2\ t . We have applied the known NLO K-factors to the signal cross sections. 
The continuum diphoton invariant mass distribution is shown as a differential cross section 
histogram with 4 GeV bins. It includes the direct photon, 1 -fragmentation and 2-fragmentation 
contributions at NLO. All samples include kinematic cuts as described in the text, and an ID 
efficiency factor e 7 = 0.8 for each photon. 
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as the case $ = 1, 1, 1, 0, so obviously one can search for optimistic cases. However, since the 
Littlest model does not define the nor even allow for educated guesses, we cannot speculate 
further. 

The situation in the [i model is somewhat more positive. Again in Fig. 4 we show the signal 
as absolute cross sections, for several choices of F\ and F 2 . The ggrj coupling is comparatively 
much larger than in the Littlest model, first because there is no 1/ \/5 hypercharge embedding 
factor penalty, and second there is the enhancement from the additional heavy states running 
in the loop. Thus, for the Golden Point, the cross section is more than an order of magnitude 
larger than the cases shown for the Littlest Higgs model. Interestingly, as a function of m v 
the product <t t; xBR(77) is nearly flat. This is due to the coincidental cancelling of dropoff 
in production cross section as m v grows, with the rise in partial width of i] — > 77. As F\ is 
taken to be closer to F 2 , the cross section becomes considerably smaller, and for F 1 = F 2 it is 
more than two orders of magnitude smaller. This is because the tirj coupling goes to zero as 
F\ — > F 2 , which not only reduces the production cross section but drives BR(77) to extremely 
small values. 

For the Golden Point at LHC, 300 fb _1 would give a 7a signal for the highest-mass point 
allowed, and the signal drops below 5a for about m v = 240 GeV. The SLHC reach would 
extend down to about m ri = 130 GeV. For the case Fi :2 = 0.6, 3 TeV, the signal is barely 5a 
for rrir] = 300 GeV, while at the SLHC one could discover the 77 in this channel down to about 
m v = 160 GeV. At the upper end of the m v ranges, the signal-to-background ratio is a decent 
S : B ~ 1/10, while at the lower end of the accessibility range it drops to S : B ~ 1/50, 
about what the SM H — > 77 inclusive search would experience. If F 12 = 1,2 TeV, the signal 
degrades further, with no observability at LHC. The SLHC, however, could have limited access, 
for about m T] > 320 GeV. The F x = F 2 cases would always be invisible. 

While we do not work out the Simple Group model case, we can expect that the signal 
would in general be a factor of 2 larger than in the fi model, based on counting the number 
of heavy quarks that run in the loop. However, the coefficients of the tii],QQi] couplings may 
have non-trivial coefficients, so this is not a rigorous prediction. 

A final note is that while we have included the NLO 1- and 2-fragmentation contributions 
to the diphoton continuum background, we have not included other sources of fake photons 
from jets, which typically are a 20 — 40% effect, depending on ultimately how well ATLAS 
and CMS can reject these events. We have also not included the NNLO signal contributions, 
which are another 20% or so. Since we further apply only K-factors for the signal with simple 
kinematic cuts, details such as 77 recoil are left out. Yet our results should still be regarded as 
conservative and optimistic since we use 4 GeV bins, which are slightly larger than necessary 
for diphoton final states. Our results are therefore a reasonable estimate of the reach of the 
LHC in these Little Higgs models with pseudo-axions at the weak scale. 

4.3 Pseudo-axion detection at a Linear Collider 

At an e + e~ collider, some models will allow for the possibility of Z* — > Hi], but we do leave 
this investigation for a later publication [31]. Here we estimate the possibility of producing 77 in 
tirj associated production. Unfortunately, even for the analogous tiH channel, the cross section 
is at most a few fb. For i], the coupling to top quarks g ttr] is proportional to v/F, so the rate 
is (v /F) 2 -suppressed. On the other hand, the rate is enhanced if m v is significantly below itlh- 
Fig. 5 shows the total cross section for this process for three choices of g t t v . 

In the mass range below 2m t the r] decays with significant BR into a bb pair. The most 
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Figure 5: The ttr] cross section at a 1 TeV linear collider, for three different values of g t t v - 
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Figure 6: Invariant bb mass distribution for e + e~ — > itbb, showing the expected number of events 
in 2 GeV bins for 1 atT 1 of data. Medium gray is the QCD contribution, dark gray is EW with 
m H = 115 GeV. The light-gray spikes are i] signals for three values of m v with g ttrj = 0.2. We 
do not include finite detector resolution effects, so e.g. the m v = 150 GeV signal is marginal. 



important background is ttbb production. In Fig. 6 we show the invariant bb mass distribution 
for signal and background, for three choices of and fixed coupling to top quarks. This plot, 
which was calculated using the programs of Refs. [35,36] shows the difficulties in detecting this 
final state: for low rj masses, the signal cross section is sizable but the QCD background is rather 
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Figure 7: Total cross section for resonant rj production at a future photon collider, with different 
anomaly coefficients as indicated by the number of non-SM particles with masses ~ 1 TeV in 
the loop. For comparison, Higgs production cross sections are shown for the SM as well as SM 
plus a heavy T quark. 

large. If the 77 mass is close to the Z or Higgs masses, the EW background is oberwhelming. For 
larger values which are well separated from the Z and Higgs, the production cross section 
rapidly decreases. In any case, the experimental analysis of the ttbb final state is nontrivial [37], 
and the achievable resolution in the invariant bb mass (with correct identification of the tt pair) 
determines the sensitivity for detecting a pseudo-axion in this final state. 

In principle, the process e + e~ — > 7* — > 777 is possible with the help of the anomaly coupling, 
but it turns out that the 667 background is generally too large for this to be viable. 

4.4 Pseudo-axions at a photon collider 

Especially for Higgs precision measurements, a high energy photon collider is expected to be 
operated at a future linear collider, by Compton backscattering laser photons off the elec- 
tron beam. At such a machine the Higgs boson(s) as well as pseudoscalar Higgs bosons can 
be produced as s-channel resonances [38]. Pseudo-axions could be produced the same way. 
The effective cross sections for resonant pseudo-axion production with different choices of the 
anomaly factors (i.e., different numbers of particles in the loop and/or different coupling fac- 
tors) are shown in Fig. 7. Since all strongly-coupled particles running in the loop are much 
heavier than the pseudo-axion itself, the cross section remains constant over a wide range of 
rrijj, while the Higgs boson shows the well-known maximum at the W threshold as well as the 
destructive interference between gauge boson and fermion loops for masses around 650 GeV. 

The decay mode 77 — > bb is typically large (although not necessarily dominant), and the 77 
manifests itself as a sharp spike in the bb invariant mass spectrum. Fig. 8 shows the bb invariant 
mass spectrum at a 200 GeV photon collider, which is optimal for studies of a light Higgs boson. 
Experimentally, the discovery is challenging since the anomaly factor must not be too small 
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and the scale F should not be too high in order for the 77 signal to be clearly visible. For larger 
r\ masses the situation is almost the same as for the linear collider: the background falls rapidly 
with increasing m b i, but the signal cross section falls at a similar rate. Fig. 9 shows the peaks 
in the m b i spectrum for larger 77 masses. For very heavy pseudo-axions, m v ~ 350 GeV, the tt 
channel opens up and one can look for peaks in this spectrum. Figs. 8 and 9 were produced 
with the programs of Refs. [39,40]. 

Muhlleitner et al. [38] studied the case of a pseudoscalar Higgs boson in a certain region of 
MSSM parameter space, in which the Higgs pseudoscalar A has exactly the same total width 
and BR(feft) as the fi model 77 at the Golden Point. They considered finite bb resolution and 
other effects, such as detector smearing of the signal, and found a significant pseudoscalar signal 
at the photon collider. 

5 Conclusions 

If a Little Higgs scenario is realized in nature, Higgs bosons at the EW scale are typically 
accompanied by new gauge-singlet pseudoscalar particles 77 that are associated with the extra 
spontaneously broken U(l) symmetry groups. If these abelian subgroups are gauged, the par- 
ticles are absorbed as the longitudinal components of extra Z' bosons. At future colliders, such 
(heavy) vector resonances can be detected by standard methods. Their indirect effects on the 
existing EW precision data already constrain the parameter space of Little Higgs models. 

Therefore, we have considered the alternative case that at least one U(l) group is ungauged, 
so that the associated NGB is physical. To avoid the limits on light axions, we require the 
presence of explicit U(l) symmetry breaking terms which give the would-be axion a mass in 
the EW range. In particular models, some of the pseudo-axion masses are calculable, while in 
other models (such as the Littlest Higgs model) is undetermined. 

Detecting these particles would be an important test of the Little Higgs model structure. At 
the LHC, one can perhaps search for them in decays of the heavy top quark partners T — > trj 
and subsequent decay 77 — > bb, or possibly 77 — > gg for extreme parameter choices or at an 
upgraded LHC or a future, higher-energy hadron collider. In some models, the decay 77 — > ZH 
is open, giving rise to extremely distinctive but probably small signals. In either case, T — > trj 
decay channels are extremely complicated and their exact utility will have to await further 
detailed work. Our work has been to point out that the presence of pseudo-axions generally 
alters the phenomenology of T quarks, often significantly and in some cases to the point of 
domination. 

Our much more interesting result is the prospect of searching for pseudo-axions in direct 
production at hadron colliders, gg — > rj, with subsequent decay to photon pairs, in the spirit 
of Higgs boson searches. While observation does not appear to generally be possible in the 
Littlest model without fine-tuned choices of the fermion U(l) v charges, in the /i model the 77 
would be visible at the LHC for larger values of m^, and over significant parameter space at the 
luminosity-upgraded SLHC. The general feature of direct production is that the more heavy 
quark partners of the model, the stronger the production rate; e.g. one would expect even larger 
rates in the Simple Group model. Unlike pseudoscalars in supersymmetry, the BR(77) is not 
suppressed by fermion couplings enhanced by tan j3 2 dominating the total decay width (tan /3 
is typically restricted to large values in supersymmetry scenarios). Instead, in general the bbrj 
couplings are suppressed relative to Higgs-like couplings, resulting in an often sub-dominant 
BR relative to the gg final state, and a non-trivial BR(77) even up to the ZH or tt thresholds. 
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The production channels at a Linear Collider (e + e~ — > ttrj) or at a photon collider (77 — > 77) 
are also promising, and we have presented results for a few cases. While the backgrounds to 
these processes appear manageable, the expected signal rates are low, so high luminosity and 
a sophisticated experimental analysis will be necessary to confirm the presence of a low-mass 
pseudo-axion. 

The pseudoscalar's nature would easily be distinguished from a standard Higgs boson by the 
absence of the WW and ZZ fusion channels, which are well-known to work over an extremely 
large mass range [41], as well as "standard" decay channels such as gg — > rj — > W + W~, ZZ [42, 
32,33]. 

To distinguish rj from a pseudoscalar Higgs in more conventional two-doublet models or 
other extended Higgs sectors, we would have to identify it as a gauge singlet. We anticipate 
that observing all possible T decay modes would aid this, but it probably would require an 
e + e~ collider which can cover the parameter space as well as prove the absence of charged and 
additional CP-even neutral partners. 
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A Heavy T quark decay partial widths 

In the Littlest Higgs model, the T partial decay widths T — > tH and T — > trj are given by 

r = l^/fo, xh) [(1 + A - 4) («i + 4) + ±K L K R x t ] (50) 

with the usual definitions x^ = mi/M? and the function 

f(xi,Xj) = \J{1~ {xi + Xj) 2 ){l - (xi - xj) 2 ), (51) 

where for H and rj the left- and right-handed couplings are 

M snf v \ .M sm 



<=0 - «£ = --- (52) 

\b / c v 

kI = i™(3' 4 = iO (£) (53) 

In the limit v/F — > where the masses of H, Z, W, and 77 can be neglected compared to Mr, 
SU (2) symmetry relates the partial decay widths into (longitudinally polarized) vector bosons 
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to the partial decay widths into a Higgs boson, and the BRs simplify to 

r(T "* tH) = 32^? = r(T tZ ° ] = 2 r(T "* bW+) (54) 

= (") 
For the \i model, the partial widths of the heavy T quark (neglecting the b quark mass) are: 

^f(*t,x h ) [(1 + af -^)(«« + 4 2 )+4«*4x t ] (56) 

^f(*t,x„) + + «J 2 )+4«2«M (57) 
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with 



A(a;,y,z) = (a; - (y + zf )(x - (y - zf ) . (61) 
Yukawa-type couplings are defined as 

t 7 "(cv - 7 5 ca)^, ^7^(1 - 7 5 )^ • (62) 

B Loop-induced couplings 

Little Higgs pseudo-axions couple to vector bosons at one-loop order via the usual triangle 
graphs in Higgs phenomenology. All fermions which get their mass by U(l) v breaking run 
in the loop. As long as these fermions are heavy compared to the axion, the loop value is 
independent of the heavy-fermion mass. It depends only on the triangle anomaly coefficient, 
i.e., the magnitudes of the effective 7777, r]W + W~ , r/ZZ and r/gg vertices, given by the mixed 
anomalies of the U(l) v symmetry with the electromagnetic, EW, and QCD gauge groups, 
respectively. We write the anomaly coefficients Cj as parameters: 

'■ (y s ~ 1 Qi ~ 

F * n F8 f . (63) 

The dual field strengths are normalized by F^ v = ^e^ upa F pa The values of the anomaly coeffi- 
cients are model-dependent and given by 

c 7 = Y. N cQ)\ F U Am )H)\ 9 -Y L ' ( 64a ) 

c 9 = EjiWiip ( 64b ) 
<? q 
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where the mass of the fermion in the loop is parameterized by m/ = XfF (for heavy particles 
A/ and g v ff are both order one, while for SM particles like the top both are order v/F) and 
F^, 2 is the function defined in [43] . We are not interested in the operators containing the W or 
the Z here because the 77 already has small production cross sections, and further small BRs 
to observable (lepton) final states makes these decays moot. The relevant couplings are given 
in Eq. (11) for the Littlest Higgs model and in Eq. (23) for the fi model. Unlike for a Higgs 
scalar, the W loop is unimportant, since the 77 does not couple to weak bosons at tree level. 
Since the W loop would otherwise contribute with opposite sign to fermion loops, there is no 
destructive interference as in the SM. 

In contrast, in Little Higgs models the masses of all new heavy fermions are expected to 
be non-invariant under the U(l) v symmetries (to make the model natural), hence they all 
contribute to the axion-gauge boson interactions with full strength. In particular, the Simple 
Group models discussed above predict heavy partners for all SM fermions. Leaving aside the 
detailed coupling structure, for the rjgg coupling the heavy quark triangle diagram value is 
approximately multiplied by N = 3 (the number of quarks) in the fi model, or N = 6 in the 
Simple Group model. For the coupling to EW gauge bosons, the heavy partners of leptons and 
neutrinos have also to be included in the loop, which we ignore here. Note also that in extended 
scalar sectors in models like the [S77(4)/S77(3)] 4 simple group, there can be enhancement effects 
by the tangent of a mixing angle analogous to that of the MSSM. The loop integrals for the 
triangle graphs are found in [43]. 

In the Littlest Higgs model, the situation is less certain since we do not know how many 
heavy fermions actually exist. Furthermore, we need the absolute U charges of the fermions, 
while the previously introduced (3 coefficients are merely differences of U(l) v charges. Here, 
we cannot predict the anomaly coefficients but have to leave them as free parameters. It is 
even possible for them to cancel for randomly chosen integer values of (3q, f3\, $2- This is largely 
irrelevant, however, as we note that the normalization factor l/\/5~ squared highly suppresses 
the anomalous partial widths. This is peculiar to the Littlest Higgs model and is due to the 
hypercharge embedding. However, judicious choice of the j3 coefficients can compensate for 
this. 

Nevertheless, it is not unreasonable to expect that the factor v/F, which suppresses the 
i] couplings compared to the corresponding CP-even or CP-odd Higgs couplings in the SM or 
MSSM, is compensated by the large weight of the heavy-fermion sector in any Little Higgs 
model. 
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